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Results are presented for the optical Hall conductivity of a Floquet topological insulator (FTI) 
for an ideal closed quantum system, as well as an open system in a nonequilibrium steady-state with 
a reservoir. The steady-state, even for the open system, is strongly dependent on the topological 
phase of the FTI, with certain phases showing a remarkable near-cancellation from pockets of Berry- 
curvature of opposite signs, leading to a suppressed low-frequency Hall conductivity, that also shows 
an anomalous temperature dependence, by increasing as the temperature of the reservoir is increased. 
Such a behavior is in complete contrast to heating, and arises because of a strong modification of 
the effective system-reservoir coupling by the laser. The Berry curvature of the Floquet modes 
is time-dependent, and its frequency components are found to control the main features of the 
high-frequency Hall conductivity. 

PACS numbers: 73.43.-f, 05.70.Ln, 03.65.Vf, 72.80.Vp 


I. INTRODUCTION 


Topological systems are characterized by a bulk-edge 
correspondence where geometric properties of the bulk 
band-structure have a precise connection with the na¬ 
ture of excitations at the edge when the system is placed 
in a confined geometry. For an integer quantum Hall 
system for example, bulk bands have a non-zero Chern 
number C, which also equals the number of chiral edge 
modeJi^. This correspondence implies that many bulk 
measurements are indirect probes of edge excitations as 
well. For example, the dc Hall conductivity at zero 
temperature, (which is equal to the Hall conductance 
measured in a multi-terminal measuremeniP) is univer¬ 
sal = axy{oj = 0) = Ce^/h, and proportional to 
the number of edge-states. Moreover, the optical Hall 
conductivity a^y (w 7^ 0), while non-universal, neverthe¬ 
less shows signatures of quantum Hall plateaus as the 
external magnetic field is variecP. In fact an all optical 
measurement such as Faraday rotation ©i;- of linearly po¬ 
larized light of frequency oj is related to the optical Hall 
conductivity as Qf{ui) ^ aCnCxyiuj) [a being the fine- 
structure constant, and c„ being a material dependent 
parameter such as the refractive index^, and has been 
used as an alternate probe of quantum Hall physics, both 
in semiconductor heterostructure^, and graphend^. 

Chern insulators are topological insulators (TIs) which 
show quantum Hall physics in the absence of a mag¬ 
netic field, where time-reversal symmetry is broken by 
introducing complex hopping amplitude^!. This can 
be achieved by the application of a circularly polar¬ 
ized lasepiiHiS, where TIs arising out of such time- 
periodic perturbations are referred to as Floquet TIs 
(FTIs(p21. The field of FTIs has grown in recent years 
because of several experimental realizations ranging from 
periodic ally s haken lattices of cold-atomic gasea^^, to 
graphen d^^b9 | and Dirac fermions on the surface of 3D 
TIs^Sl under external irradiation, and also arrays of 
twisted photonic waveguided^. In fact FTIs are ex¬ 
tremely rich, showing different topological phases as the 



FIG. 1. Optical Hall conductivity (main panel) with tem¬ 
perature dependence of the dc Hall conductivity (inset) for 
laser frequency Q. — 5.0th- The laser amplitude and Chern 
number are upper-panel: Aqo = 1.5, C = 1, lower-panel: 
Aon = 0.5, (7 = 3 . 


amplitude, frequency, and polarization of the periodic 
drive is variecP^H^. 

In this work we present results for the optical Hall con¬ 
ductivity of FTIs by accounting for the fact that these 
systems are far out of equilibrium. To this end, the pre¬ 
cise relaxation mechanisms need to be specified as they 
sensitively affect results. Thus we present results for two 
rather different cases, one where the system is closed and 
the laser is switched on as a quench, while the second 
is when the system is coupled to an ideal reservoir, and 
so the steady-state loses memory of its initial state. Yet 
being a driven dissipative system, the steady-state is far 
from thermal, resulting in some unusual behavior for the 
Hall conductivity, which is sensitive to the topological 
phase, and cannot be interpreted as simply “heating”. 
In fact we show that the laser parameters, and in partic¬ 
ular the topological phase, strongly affects the effective 
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FIG. 2. Quasi-energies in a cylindrical geometry highlight¬ 
ing edge states. The laser frequency is = 5.0th, while 
the laser amplitude and Chern number are (upper panel) 
Aoa = 1.5, C = 1 and (lower panel) ^oa = 0.5, G = 3. The 
latter shows edge-states both at the center and the edges of 
the FEZ. 
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FIG. 4. (color online) Upper panel: Time-averaged Berry 
curvature over 4 BZs for fl = 5.0th, A^a = 0.5, C = 3. Lower 
panel: Gontour plot for the steady-state population difference 
for reservoir temperature T — 0.01th- Note the almost com¬ 
plete population inversion between regions around the Dirac 
points and the circles. 



FIG. 3. (color online) Time-averaged Berry curvature over 4 
BZs for Q, = 5th, Aoa = 1.5 corresponding to a Ghern number 
of G = 1. The structure is rather smooth, moreover in the 
Fourier decomposition of Fkd{t), two harmonics (F’™^'^’^) are 
dominant. 


system-reservoir coupling, and therefore the steady-state, 
and its Hall response. 

The advantage of the optical Hall conductivity is that 
it can be measured both in traditional multi-terminal 
measurement^2lHlH!^ ag -^vell as completely optically via 
the Faraday rotation, the latter being more suitable for 
present day experiments. In fact even the dc limit can 
be studied without leads, in optical lattices, by monitor¬ 
ing the transverse drift in time-of-flight measurementJ^^. 
Moreover while it is the dc (w —?► 0) limit of the Hall con¬ 
ductivity, combined with an ideal electron distribution 
function, that is related to a topological q uantit y, namely 
the Chern number of the Floquet band^^HHl^ -^^e show 
that even at non-zero frequencies, some features of the 
topological nature of the system survive. In particular, 
the Berry-curvature of Floquet bands is time-dependent, 
and we show that various frequency components of this 
Berry-curvature can be directly probed in the optical Hall 
conductivity. 

The paper is organized as follows. The model is in¬ 
troduced in Section [n| and the derivation of the Floquet- 
Master equation for the open system is outlined (with 


the optical Hall conductivity is derived using a linear- 
response approach, while in Section |IV| results for the 


some details relegated to Appendix In Section hlh 
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optical Hall conductivity are presented. Finally in Sec¬ 
tion]^ we present our conclusions. 


II. MODEL 


Our model is graphene irradiated by a circularly po¬ 
larized laser, and also coupled to a phonon bath. The 
Hamiltonian is, H = He\+Hp\^-\-Hc where (setting h= 1) 
Hei is the electronic part, 


^^el= —th ^ 


„t J 

^kB 


)( 


0 






CkA 


0 J \CkB 


( 1 ) 


where = X)i=i 2 3 Si are the nearest- 

neighbor unit-vectors of graphene, a is the lattice spac¬ 
ing, and A is the circularly polarized laser of amplitude 
Ho and frequency H, A^it) = 9{t)AQ cos fit; Ay (t) = 
—0(t)Ho sin Hi, which we assume has been suddenly 
switched on at time t = 0. We consider dissipation aris¬ 
ing due to 2D phonons Ffph = J2q,i=x,y 
are coupled to the electrons as 




which 


Hr = 




/ ^ ’ ^tyiy'Okrj' ( 2 ) 

kqa,(7' —A,B 


tj is a pseudospin label denoting the A, B sub-lattices, 
and the electron-phonon coupling Mg\_pi^{k,q) = 


A. 


',kq 


bx,q ) ; ^y,kq 


is off- 
While we will 


diagonal in pseudo-spin spac^2SIlI]_ 
consider only inter-(quasi)-energy-band transitions 
driven by long-wavelength optical phonons, our re¬ 
sult showing that the effective matrix-elements of the 
electron-phonon coupling is modified by the laser, lead¬ 
ing to nonequilibrium steady-states that depend strongly 
on the topological phase of the FTI, is rather general, 
and will carry over to other kinds of system-reservoir 
couplings. 

We assume that initially the system is in the ground 
state of graphene |'l'(t = 0)) = Hs; IV'in.k), and we time- 
evolve the system according to H. For the closed system 


(He = 0), this time-evolution can be studied exactly, 
while in the presence of phonons, we solve the problem 
by assuming a weak electron-phonon coupling, and a very 
fast bath s o that a Floquet-Markov approximation can be 

-EH. 


Let us denote W{t) to be the full density ma¬ 
trix obeying dW{t)/dt = —i[H,W{t)], which in 

the interaction representation becomes, Wi{t) = 

where = 

rife Uk{t, t') is the time-evolution operator for the 
electrons under a periodic drive, and Uk{t,to) = 
Ea=«.de"*"'=“^‘"*“^|<('fea(i))(<('fe-a(io)| with |(^fe.a(t)),efe„ 
being the quasi-modes and quasi-energies that obey 


iLei - idt 


l^fea) 


^ka\4>ka)- We restrict the quasi¬ 


energies to a Floquet Brillouin zone (FEZ), —ft/2 < 
Cka < fl/‘^i where the two quasi-energy levels are la¬ 
beled €ku > Cfed- The electron reduced density matrix, 
obtained from tracing over the phonons, = TrphVF, 
has the following form in the interaction representa¬ 
tion, Wf{t) = Pk,afi{tMka{t)){(l)kl3{t)\. For a 

quench switch-on, the physically relevant initial condi¬ 
tion is given by the overlap between the initial state and 
the Floquet modes at t = 0, pk,ai 3 {t = 0) = = 

(^fect(O)|'!/'in,fe)('0in,fe|^^'fe/3(O)). Under assumptions that the 
bath is Markovian, always stays in thermal equilibrium 
at temperature T, and that the electron-reservoir cou¬ 
plings small in comparison to all quasi-energy level spac- 
ings^S, the density matrix obeys the equation of motion, 
Pfe.aa(t) = -E/3=«.d^L;/3/3Pfc:/3/3(0, where are 

the scattering rates. The steady-state solution is: 


Pk,oia{^ — — Pk 

\lL. 

Pku — 

\^uu.dd\ + \^UX 


dd\ 


, Pkd = 1 — Pku 


( 3 ) 

( 4 ) 


where the scattering rates for a constant electron-phonon 
coupling (\x,qk = \,qk = A) with a broad phonon den¬ 
sity of states Dph are. 


rfc 

^UU^UU 




^{~^kd ~\~ ^ku + 'klQ.)N{ — Ckd + ^ku + bQ) -h O^Ckd ~ ^ku ~ ‘kl^)N{ekd ~ ^ku ~ Tl^) 


rk 

^uu,dd 


E "J 


n=int 


^i^^kd ^ku 


nfl)N{ekd - e.ku - nfi) + 9{-ekd + Cku + nft)N{-ekd + efe« + nfl) 


( 5 ) 

( 6 ) 


Above 9{x) is the step-function, N{x) = l/(e^/^ - 1) {(j)ka{t)\c\^Cks\(/kp{t)) = Ln^int e“^*™afe,ap¬ 

is the Bose function, N = 1 -|- A^. The matrix el¬ 
ements for scattering of electrons between the quasi¬ 
energy levels Ckd £ku + nfl by phonon absorption or 

emission are MJ} = {'m-//k,ud^dk,du) > ^^lere 
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III. DERIVATION OF THE OPTICAL HALL 
CONDUCTVITY 


We will now explore the optical Hall conductivity a^y 
of the open and closed system, where for the latter 
Pka = is & response to a weak perturbation 

applied over and abo ve the laser, and is computed using 
linear-response theorj^HHl particular the current in 
the direction z, in response to a weak electric field E{t), 
is ji{t) = f dt'aij{t,t')Ej{t'). In general cjij{t,t') is not 
time translationally invariant due to the time-periodic 
perturbation. We present results after time-averaging 
t + t' over a laser cycle Tq_ = 27r/fl, and then Fourier 
transforming with respect to the time difference t — t'. 
Before presenting explicit expressions for the <Jxy, note 
that due to the time-dependence of the Floquet modes, 
the Berry curvature 


Fkd{t) = 2Im 


{dy(j)kd{t)\dx(j)kd{t)) 


( 7 ) 


is time-dependent. Yet its integral over the BZ is time- 
independent and gives the Chern number 

C=^ [ d^kEkdit) ( 8 ) 

Jbz 

Defining the Fourier transform of the Berry “vector po¬ 
tential” , 

1 r) 

a natural object that appears in the optical Hall conduc¬ 
tivity is. 


Fl^=^ 


-^uxd-^dyu ■^uyd'^dxu 


( 10 ) 


which can be thought of as a Fourier decomposition 
of Fkd{t) and is such that X)m=int the Berry- 

curvature time-averaged over one cycle of the laser, 


E 

m=int 


2 

% 


(•Tn 


dtlm. 


{dy(l)kd{t)\dx(j)kd{t)) 


( 11 ) 


We find that the optical Hall conductivity at steady-state 
may be written as. 


(Jxy{uj)= E 

m—int 

where cr™ depend on FJ^ as (i5 = 0+), 


( 12 ) 




^ku ^kd TTtFl 


f: 


- {f-ku - <^kd - rnilf - 2i6u! 


- (efe„ - ekd - mny 




Pkd Pku 


(13) 


In the low frequency (dc) limit, this reduces tcP^, 

cra;y(w = 0) = —— / dykFkd[pkd- Pku] (14) 
27r/i Jbz 

The optical Hall conductivity depends on the occupa¬ 
tion probabilities Pka=u,d and we refer to the “ideal” 
limit as one where only one Floquet band is fully oc¬ 
cupied \pkd — Pku\ = 1, so that the dc Hall conductivity 
is = 0) = Ce^ jh. We will show below that while 

controls the low frequency Hall conductivity, the 
individual FJJ control the high-frequency Hall conductiv¬ 
ity which show enhancement (side-bands) in the vicinity 
of w ^ \eku — ^kd — rnn] where k are those points in 
the BZ where F"™ are peaked. We will present results for 
Re [axy{oj)], and choose a disorder broadening 5 = H/IOO. 


IV. RESULTS 


The system shows a series of topological phase transi¬ 
tions as the laser amplitude or frequency is varied, with 
the topological transitions involving level crossings at the 
center and/or boundaries of the FB^^. Moreover while 
topological phases with non-zero Chern number are pos¬ 
sible both for laser frequencies off-resonant (H > Qth) and 
resonant (H < Qth) with the electronic states, the topo¬ 
logical phases may be rather different for these two cases, 
as for the latter, the laser can create an effective band 
inversiorfi^, strongly modifying the Berry curvature. We 
find that the inelastic matrix-elements MJ} also depend 
on the laser parameters, by being strongly peaked around 
n = riQ, where no depends on the resonance condition 
(i.e., whether the laser is off-resonant, or single-photon 
or two-photon etc processes dominate). 

In particular, in the limit of large (off-resonant) laser 
frequency, and small laser amplitude {Agath/^ 1), 
the term dominates, and the resultant distri¬ 

bution is very similar to the conventional one pku 


1 / 


exp {eku - ekd)/T + I 


and thus for this case the 


reservoir ’’cools” the system, giving a Hall response that 
approaches the ideal limit as the reservoir temperature 
is lowered^. In contrast, as discussed below, topologi¬ 
cal phases where the laser frequency is smaller than the 
band-width can give rise to dominant matrix elements 
where ng are strongly quasi-momentum depen¬ 
dent, resulting in unusual Hall response. 

Fig. n shows axy along with its temperature depen¬ 
dence, lor laser frequency H = 5.0t/i and two different 
amplitudes Ago = 1.5, 0.5, corresponding to two differ¬ 
ent topological phases (7 = 1,3 respectively. The quasi¬ 
energy spectra for these two phases in a cylindrical ge¬ 
ometry highlighting edge states is shown in Fig. while 
the Berry-curvature of the phase with Ago = 1.5, (7=1 
is shown in Fig. and that of the phase with Aga = 
0.5, (7 = 3 is shown in Fig. 

One finds that the Chern insulator for Ago = 1.5, (7 = 
1 is a conventional one (ng = 0), with the bath effectively 
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cooling the system in that the Hall conductivity increases 
as the temperature of the reservoir is lowered. In contrast 
the case of A^a = 0.5, C = 3 is quite different because 
the low temperature Hall conductivity is small, almost 
zero, and has an anomalous temperature dependence in 
that it actually increases as the reservoir temperature is 
increased. In fact the entire region of H = 5t/j,0.01 < 
A^a < IjC = 3 constitutes the same topological phase, 
and shares this behavior of low dc Hall conductivity, and 
non-monotonic temperature dependence. 

Fig. a shows the time-averaged Berry curvature and 
the population imbalance for the above phase {A^a = 
0.5, (7 = 3). Besides the characteristic peaks at the Dirac 
points that one expects in Chern insulators described by 
the Haldane model, Fkd for these laser parameters also 
shows sharp circular rings. While the Dirac points con¬ 
tribute to a Chern number of 1/2, the circular rings give 
a Chern number of 2, so that in total (7 = 3 (since the 
BZ contains two Dirac points and a ring). The steady- 
state distribution function for this case shows that there 
is effectively a population inversion between the regions 
around the Dirac points where the electrons are primarily 
in the ’’down” level, and the regions within the circles, 
where the electrons are primarily in the “up” level. This 
population inversion arises due to the structure of the 
the matrix elements Tf/*” where in the vicinity of the 
Dirac points is dominant, while within the circles 
the ~ 0, while are dominant. This im¬ 

plies that at low temperatures where phonon absorption 
is suppressed, the system can only relax via phonon emis¬ 
sion from tkd -l- D —> Eku within the circles, whereas in 
the regions around the Dirac cones, the phonon emission 
processes cause relaxation from Cku —>■ Cfed- Since the Hall 
conductivity involves integrating over the entire BZ, this 
results in a low Hall conductivity. 

One way to differentiate between this subtle matrix- 
element effect, and simply heating, where the latter will 
also give a low Oxy, is by studying the temperature de¬ 
pendence. As the temperature is increased, one finds that 
the Hall conductivity actually increases (opposite to the 
heating case). This happens because the quasi-energy 
level separation near the Dirac points and around the 
circles is not the same, with the former being larger than 
the latter. Raising the temperature excites electrons to 
the ’’down” level near the circles, while excitations to the 
”up” level near the Dirac cones do not occur until higher 
temperatures. This results in a low temperature region 
where the Hall conductivity increases with temperature 
because in this regime as the temperature increases, the 
effective population imbalance between the two Floquet 
bands with opposite Berry curvature increases. This 
qualitative behavior is not special to these parameters, 
but holds for the entire topological phase where the only 
difference is the size of the circle in Fig. and more 
generically is a signature of the fact that more than one 
scattering rate plays a dominant role in transport. 

We now discuss this strong matrix element depen¬ 
dence in another phase of the FTI. We choose a much 
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FIG. 5. Quasi-energies in a cylindrical geometry highlight¬ 
ing edge states. The laser frequency is D = 0.5th, while 
the laser amplitude and Chern number are (upper panel) 
Aoa = 10, (7 = 0 and (lower panel) Aqu = 5,(7 = 6.0. It 
is the top phase (C — 0) whose Hall response is discussed in 
detail in the text. 



FIG. 6. (color online) Time-averaged Berry curvature over 
4 BZs for D = 0.5th,Aoa = 10 corresponding to a Chern 
number of (7 = 0. 
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FIG. 7. (color online) Sign(no) in the dominant matrix ele¬ 
ment for = 0.5th, Aqu = 10, C = 0. The sign of the 
population imbalance {pku — pkd) entering the Hall response 
follows this pattern. 


lower laser frequency, but larger laser amplitude (fl = 
0.5th, A^a = 10, C = 0) where one expects several band 
inversions as resonances involving n = 1... 12-photon 
processes are in principle allowed, although how visible 
these are, depends on the laser amplitude. The corre¬ 
sponding quasi-energy spectra and Berry-curvature are 
shown in Fig. and respectively. 

For Aqq = 10, MJt has contributions from many n 
where the dominant matrix element 44^° is such that 
—4 < no < 4 in the BZ. Fig. shows how the sign of no 
varies through the BZ, where for no = 0 or a positive in¬ 
teger, the population is primarily in the ’’down” level at 
low temperatures, while no-negative gives a population 
primarily in the ”up” level. The sign of the population 
difference follows the pattern in Fig. leading to a dc 
Hall conductivity that is not only small, but also has a 
non-monotonic temperature dependence. Changing the 
laser parameters in such a way as to stay within the same 
topological phase, simply changes the size of the pock¬ 
ets of max(no) in the BZ, with no crossings or no new 
pockets appearing. Note that this particular phase is not 
a usual Chern insulator as C = 0, yet the system does 
support edge states. Moreover, while we are studying the 
system in the bulk, and edges do not enter explicitly in 
our Kubo formula approach, yet we do pick up a non-zero 
Hall response coming from the non-zero Berry-curvature 
of the bulk bands. 

While so far we have been mainly discussing the low 
frequency Hall response, we now make some general ob¬ 
servations about the high frequency response. We find 
that for small amplitudes Aqq ~ (2(1), Oxyiuj) is satu¬ 
rated by w ~ H (Fig. [^. On the other hand, for larger 
amplitudes, like the case just discussed {A^a ~ 0(10)), 





FIG. 8. axy{u!) for Q — 0.5th,A(,a = 10, C = 0. In compari¬ 
son to the case of O = 5.0th, Aoa = 1.5, 0.5, many side-bands 
are visible, and the response is sensitive to the distribution 
function over a larger frequency range ta. 


axy is non-zero over a larger range of frequencies (see 
Fig. §, with the position and magnitude of the side¬ 
bands depending on the structure of F^. This is high¬ 
lighted in Fig. 1^ where the structure in the total optical 
Hall conductivity plotted in Fig. has been decomposed 
into those arising from separate Fourier components of 
Fk- 

To understand the location of the peaks in w, it 
is convenient to look again at the low amplitude case 
(n = 5.0th, A^a = 1.5&0.5). Here only are dom¬ 

inant, but while for one (^oa = 1-5), FJA^^ is peaked 
at €ku — ^kd ‘C n, for the other {A^a = 0.5), is 

peaked at etu — t-kd ~ In fact are re¬ 

spectively the Dirac cones and circles in Fig. Thus a 
side-band appears at w ^ \eku — t-kd — fl] 0.7512 for 
the case of Aqu = 1.5, C = 1, while no such side-band is 
visible for A^a = 0.5, (7 = 3. Similar arguments can be 
utilized to understand the peaks in the high frequency 
response for the large amplitude case shown in Figs. 
and Thus while the optical Hall conductivity at non¬ 
zero frequency cfxyiuj ^ 0) is not a topological quantity, 
yet it is sensitive to the structure of the Berry curvature 
both in momentum and time. 


V. CONCLUSIONS 

In summary we have presented results for the opti¬ 
cal Hall conductivity which can be measured in an all 
optical measurement such as Faraday rotation, as well 
as using leads. The low frequency Hall conductivity, 
even in steady-state with a reservoir, is remarkably sen¬ 
sitive to the topological phase, which is communicated 
through the structure of the system-reservoir matrix el¬ 
ements, resulting in steady-states with anomalous tem¬ 
perature dependence, and also opening up the possibility 
of manipulating the steady-state by suitable reservoir en- 
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FIG. 9. a^y for = 0.51;,, ^ 0 “ = 10, G = 0 and for an ideal 
distribution pku — pkd = 1. The peaks are at cu ~ |efe„ — tkd — 
mfll where the location of k is determined by the peaks in 

Fr- 


We presented detailed results for the optical Hall re¬ 
sponse in three different phases of the FTI. One is where 
the Chern number coincided with the number of edge- 
states which were all in the center of the FEZ (H = 
5th, Aqq = 1.5, C = 1), the other where the Chern num¬ 
ber does not equal the number of edge-states but rather 
equals the difference in the number of chiral edge-states 
above and below the Floquet band (H = 5th, A^a = 
0.5, C = 3), while the third was one where (7 = 0 
(H = Q.5th,AQa = 10,(7 = 0), so that the band was 
not a conventional Chern insulator, yet the system sup¬ 
ports edge states with the number of chiral edge modes 
above and below the band equal to each other (hence 
(7 = 0). An interesting future direction of research is to 
compare the results presented here which were based on 
a bulk Kubo formula treatment, and which is sensitive to 
the bulk Berry-curvature, with a Landauer transport ap¬ 
proach through a short sample where the leads determine 
the occupation of the Floquet states. 


gineerinj2lllll^ for example by applying strain fields that 
modify the symmetries of the electron-phonon coupling. 
We also find that the high-frequency response acts as 
a spectroscopic tool for the Fourier components of the 
time-dependent Berry curvature. 
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Appendix A: Derivation of Floquet-Master equation 


We assume that initially the system is in the ground 
state of graphene |'L(t = 0)) = Hfc IV'in.k), and we time- 
evolve the system according to H. For the closed system 
{He = 0), this time-evolution can be studied exactly, 
while in the presence of phonons, we solve the problem 
by assuming a weak electron-phonon coupling, and a very 
fast bath so that a Floquet-Markov approximation can 
be made . Th is approach has been discussed in detail 
elsewher d^^ l ^^ l, however for completeness we briefly high¬ 
light the main steps. Let W{t) be the full density matrix 
obeying = —i[H,W{t)]. In the interaction rep¬ 

resentation, Wi{t) = e*'^p‘'*C/j^^(^, 0)LF(t){7ei(t, 0)e“*^p‘'‘, 
where Uei{t,t') = is the time-evolution op¬ 

erator for the electrons under a periodic drive where 
Uk{t,to) = with 

\4>ka{t)) T ^ka being the quasi-modes and quasi-energies 


that obey 


iJei - idt 


\4^ka} 


£ka\4>ka)- We will follow 


the convention of restricting the quasi-energies to a Flo- 
quet Brillouin zone (FBZ), —17/2 < eka < f^/2, and label 
the two quasi energy-levels such that Cku > ^kd- Defining 
the electron reduced density matrix as the one obtained 
from tracing over the phonons, = TrphVF, at 0{Hl) 
we need to solve, 

dwf r* 

^=-Trph / dt'[Hej{t),[He,i{t'),Wi{t')]]{Al) 
at Jto 

We assume that at the initial time tg, the electrons 
and phonons are uncoupled so that W{to) = W(f'(to) G) 
LFP^(to), and that initially the electrons are in the post- 
quench state |^'(t)) = Jlfc Dfc(t, 0)IV'in.k)- Thus, W^\t) = 
= rtfc where 

Q,^—u,d 

with = (V>fcc<(0)|V'in.fc)(V’in.fc|'Vfe/3(0)). Since the 

phonons are an ideal reservoir that stay in thermal equi¬ 
librium at temperature T at all times, we write Wi{t) = 
Wf{t) G) e“-^p''/’^/Tr [e“^p**/^] (we set ks = !)• 

The most general form of the reduced den¬ 
sity matrix for the electrons is Wf{t) = 
rifc Ea/3 Pk,a0{t)\4>ka{t)){4'kp{t)\ where in the absence of 
phonons, pk,ap = P^^p^ and are time-independent in 
the interaction representation. The last remainiM step 
is to make the following three approximationgSSnSl (a). 
Pk,ap are slowly varying as compared to the characteris¬ 
tic time scales of the reservoir (Markov approximation), 
(b) they are also slowly varying in comparison to the 
laser frequency so that the scattering rates can be aver¬ 
aged over one cycle of the laser (modified rotating wave 
approximation) and, (c) we are away from any topologi¬ 
cal transitions so that the spacing between quasi-energy 
levels is large as compared to the coupling to the bath 
{\^ku - ikd - nil] > DphA^, |nf7| > DphA^), under such 
conditions the off-diagonal matrix elements (pkud) are 
small and can be dropped. These approximations lead to 
the rate equation, pk,aa{t) = -J2p=u,d^ta;ppPk,i3p{t) 
where pp the scattering rates. 









